The onset of jamming as the sudden emergence of an infinite fc-core cluster 
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A theory is constructed to describe the zero-temperature jamming transition of repulsive soft 
spheres as the density is increased. Local mechanical stability imposes a constraint on the minimum 
number of bonds per particle; we argue that this constraint suggests an analogy to fc-core percolation. 
The latter model can be solved exactly on the Bethe lattice, and the resulting transition has a mixed 
first-order/continuous character reminiscent of the jamming transition. In particular, the exponents 
characterizing the continuous parts of both transitions appear to be the same. Finally, numerical 
simulations suggest that in finite dimensions the fc-core transition can be discontinuous with a 
nontrivial diverging correlation length. 

PACS numbers: 64.60.Ak, 64.70.Pf, 83.80.Fg 



Understanding a continuous phase transition is tan- 
tamount to determining the universality class to which 
it belongs. In contrast, understanding the nature of 
a discontinuous change of phase requires a detailed 
study of the system at hand. Under normal circum- 
stances Q, the two categories are mutually exclusive. 
However, there are a few examples of continuous tran- 
sitions that exhibit characteristics of first-order transi- 
tions HH&EEHBHII3 In this Letter ' we wil1 

present arguments that the jamming transition in sphere 
packings 0, 13 El belongs to this class and can gen- 
uinely be described as both continuous and discontinu- 
ous. Indeed, we will identify the minimal physics needed 
to capture the nature of the transition by analogy to the 
It-core percolation model, and show by exact calculation 
that the latter model has a true mixed transition of this 
type with similar exponents at the level of mean-field 
theory. We also present numerical evidence that fc-core 
models can still exhibit mixed transitions in finite dimen- 
sions. We remark that , st arting from a different vantage 
point, Toninelli, et al. |l4| have arrived at a model of the 
A;-core type and have reached similar conclusions about 
the nature of the transition in their studies of kinetically- 
constrained models. 

Numerical studies 0, 0, El of sphere packings at 
zero temperature suggest that there is a packing den- 
sity <j) c (Point J) where the onset of jamming is truly 
sharp; i.e. the static bulk and shear moduli vanish 
for 4> < 4> c and are nonzero for <f> > <f> c . This transi- 
tion exists for spheres that repel when they overlap and 
otherwise do not interact. For small <f>, particles eas- 
ily arrange themselves so as not to overlap with any 
other particle and hence the total potential energy is 
V = 0. As is increased, there is a particular value of 
cj) c above which the particles can no longer "avoid" each 
other and V becomes nonzero. The average coordination 



number (the average number of overlapping neighbors 
per particle) is Z = for <f> < </> c . As <f> approaches 
(f> c from above, however, the behavior is very different: 
< Z >« Z c + Z ((j> - <j> c f, where (3 = 0.49 ± 0.04 [l^. 
Moreover, the singular part of the shear modulus van- 
ishes with the exponent 7 = 0.48 ± 0.05 01 an d recent 
simulations by Silbert, et al. Il3l fi nd a diverging length 
scale exponent v = 0.24 ± 0.030, El 

These numerical results imply that the transition at 
Point J has characteristics of both types: certainly there 
is a discontinuity in the average coordination number, 
(Z), but as the transition is approached from the or- 
dered (jammed) phase, it exhibits the typical singulari- 
ties associated with continuous transitions; (Z) tends to 
its limiting value with a nontrivial power-law and there 
are divergent length scales. 

We will now present arguments that the Point J transi- 
tion is indeed a mixed transition, many aspects of which 
can properly be understood by analogy to a relatively 
simple model called "fc-core percolation" (sometimes also 
called "bootstrap percolation" Let us start with 

an informal discussion of the essentials of the jamming 
model. Clearly, a jammed packing of spheres at T = 
must be mechanically stable. For a sphere in d dimen- 
sions to be locally stable, it must have interactions (i. e. 
overlap) with at least d+l neighboring spheres • Ev- 
idently, spheres with fewer than d+l overlapping neigh- 
bors do not contribute to the formation of a jammed 
structure and thus are irrelevant. Thus we may envision 
the mechanics for a system below the jamming threshold 
density as its energy is lowered towards the minimum: al- 
though large clusters of overlapping particles may happen 
to form, those at the boundary of the cluster are unstable 
and will move away, further lowering the energy. This in 
turn exposes secondary particles, who are in turn forced 
to move away, and so forth until the cluster dissolves. At 
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high density the situation is more complicated. However, 
it is still true that all particles that do contribute to the 
jammed structure must have at least d + 1 overlapping 
neighbors that are not "irrelevant", and each of these 
overlapping neighbors must have at least d + 1 overlap- 
ping neighbors that are not irrelevant, and so on. In other 
words, only particles that survive this entire hierarchy of 
irrelevance can contribute to the jammed structure. 

These considerations are suggestive of the fc-core per- 
colation model, defined as follows. Consider a regular 
lattice of coordination number Z max and some integer k 
with 2 < k < Z mM . Initially, sites are independently 
occupied with probability p. In the first stage, all occu- 
pied sites with fewer than k neighboring occupied sites 
are eliminated. Then, this decimation process is applied 
to the surviving occupied sites, and so on, until all sur- 
viving sites (if any) have at least k surviving neighbors. 
Thus, at the end of this process, every surviving site has 
at least k neighbors, all of whom in turn have at least k 
neighbors, etc. The surviving sites are called the fc-core 
and phases of the model are determined by the presence 
or absence of an infinite cluster of these survivors. 

The overall analogy between the two models is self- 
evident. The initiating density p corresponds the the 
packing fraction <f>, k corresponds to d+ 1 and Z max to the 
so-called kissing number, that is the maximum number of 
equivalent hyperspheres in d dimensions that can touch 
a central one without overlaps. (In d — 2, Z max = 6, and 
in d = 3, Z max = 12.) 

In the mean-field theory of fc-core percolation, i.e. the 
Bethe lattice and infinite-range complete graph models, 
it is well-established that the order parameter undergoes 
a discontinuous jump at threshold I la. accompa- 
nied by a square-root singularity |l6l Il8l Tl9j. However, 
the fact that the latter was indicative of a critical phe- 
nomenon has heretofore been underemphasized; in par- 
ticular, a divergent length scale had not been identified. 
Below we will show that, at least on the Bethe lattice, 
there is indeed critical behavior in the sense of a diverg- 
ing susceptibility and correlation length and the various 
exponents are in (rather dramatic) agreement with their 
counterparts in the Point J simulations. 

The fc-core percolation model can be solved exactly on 
the Bethe lattice. We begin by considering the half-space 
Bethe lattice, for which we derive recursion relations for 
quantities at level n + 1 in terms of quantities at level 
n. All occupied sites at level of the half-space Bethe 
lattice belong to what we call the "deep core." We keep 
track of two quantities, the probability of belonging to 
the deep core at level n + 1, T^ l7 which requires the 
site at level n + 1 to have at least k neighbors at level 
n that belong to either the deep core or to what we call 
the "corona." To be in the corona at level n + 1, a site 
must have exactly k — 1 neighbors at level n that belong 
to the deep core or the corona. We denote by Q^+i the 
probability of belonging to the corona at level n + 1, and 



by ^n+i = Tn+i + ^n+i tn e probability of belonging to 
either the deep core or the corona. The deep core will 
necessarily be part of the fc-core; we need to keep track 
of the corona because when two half-spaces are glued 
together to form the full Bethe lattice, corona can be 
converted to fc-core. The recursion relation is 

ee P n£» r (rf ). (1) 

In the limit of large n, = = T HS . Clearly, 

pHS = o is always a solution. However, there can be a 
nontrivial solution for p exceeding some p c . For k > 3 

T HS ~a + b(p-p c ) 1/2 (2) 

At the transition, the curve pn^_i(r HS ) is just tangent 
to T HS , i. e. PcIi'^iT^) = 1. 

The average coordination number, susceptibility and 
correlation length exponents, which are needed for the 
comparison to sphere packings, must be calculated on 
the full Bethe lattice, obtained by connecting two half- 
space lattices. The resulting probability of belonging to 
the fc-core, K, is given by K = T HS + $ H Sr HS . This has 
the same behavior as T HS in Eq. 2. The average number 
of occupied neighboring sites per occupied site (i .e. the 
average coordination number) also behaves in the same 
fashion as T HS (Eq.0). It jumps from zero for p < p c to 
(Z) w Z c +Zo(p— Pc) 1 / 2 forp > p c , in excellent agreement 
with the numerical results for sphere packings [12j. 

The susceptibility is the sum of correlation functions, 
T£ jm , connecting levels I and m of the Bethe lattice, and 
has the form x — S n (^roax — l)' l/r o,n- We consider two 
different correlation functions: Tq n represents the prob- 
ability that both level and level n are connected to the 
deep core, while Tq n represents the probability that levels 
and n are connected to each other via the corona [20| . 
This latter probability can be derived by considering the 
chain of sites connecting a given site at level to a site at 
level n. To belong to the corona, each site along the chain 
must be occupied and have exactly k neighbors (including 
the two adjacent sites along the chain) connected to the 
deep core or corona. The corresponding probability for 
each site is 8 = P^T^ 2 ) (r HS ) fc - 2 (l -r HS ) z — ~ k . The 
final probability r * n therefore scales as 9™. When this 
is summed over n, it yields a susceptibility exponent of 
7* = 1/2. By somewhat more complicated reasoning, the 
dominant contribution to T^ n scales as n0™ ; this leads to 
7# = 1. Note that r * n measures the size of the corona, 
which is the region that can be converted into fc-core or 
not, depending on the state of only one site. This is the 
source of cooperativity underlying the transition. 

Another way to compute a susceptibility is to calculate 
the response to a perturbation; in the case of percolation, 
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this corresponds to the addition of low-probability "short 
routes to the infinite cluster." In our system we have 
done this in two ways: first by providing a small fraction 
of additional random sites with fc — 2 occupied neighbors 
that are connected directly to the deep core plus corona, 
and secondly by declaring a small fraction of sites to be in 
the deep core plus corona regardless of their connectivity. 
Both prescriptions yield 7* = 1/2, as well as a "magnetic 
field" exponent of 6* =2. 

The correlation length corresponding to both suscep- 
tibility exponents diverges with the exponent v* = 1/4 
(although there is a logarthmic difference between the 
two). We use the usual embedding of a Bethe lattice in 
Euclidean space HU to arrive at this result. One would 
expect the usual mean field relation v = 7/2 to hold; the 
exponent 7^ = 1 may be an artifact of the Bethe lattice. 
However, we also obtain v# = 1/2 by looking at how 
quickly the order parameter approaches its bulk value as 
a function of distance from the boundary. 

The exponents /3 = 1/2, 7* = 1/2, and v* = 1/4 
are in excellent agreement with numerical simulations of 
particle packings near Point J. However, these simula- 
tions are carried out in 2 and 3 dimensions while the fc- 
core calculations correspond to infinite dimensions (the 
mean- field limit). This raises the question of whether 
the mixed nature of the fc-core transition can survive in 
finite dimensions. Most studies have focused on particu- 
larly simple systems such as the 2d square and triangular 
lattices J22j, some 3d cubic lattices |23j and hypercubic 
lattices |9(- For these simple systems, the transition falls 
into one of two categories: Either the transition is con- 
tinuous or it does not occur until p — 1. Systems that 
exhibit continuous transitions all contain self-sustaining 
clusters, i. e. clusters that are finite and yet survive the 
decimation process. For example, for k = 3 on the 2d 
triangular lattice, the smallest self-sustaining cluster is a 
fully-occupied hexagon and the fc-core transition appears 
to correspond to ordinary percolation of these hexagons 
[22l | . Systems that fail to exhibit a transition below p = 1 
apparently contain "unstable voids" 0, 124I Ejf that lead 
to decimation of the entire population whenever p < 1. 

We regard percolation of self-sustaining clusters and 
unstable voids as "artifacts" of simple fc-core models. In- 
deed, for jammed sphere packings these two effects can- 
not arise. First, self-sustaining clusters of overlapping 
particles are forbidden due to the repulsive nature of 
the interactions between particles. Likewise, voids (i.e. 
collections of floaters, or particles with no overlapping 
neighbors) cannot grow because of the force constraints 
and because floaters can shift around, but cannot actu- 
ally disappear. To see this, consider the interface be- 
tween a void and the surrounding sea of particles with 
at least fc overlapping neighbors. For a sufficiently large 
void, a particle on the boundary with at least fc over- 
lapping neighbors will inevitably experience a nonzero 
net force into the void, since floaters provide no compen- 
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FIG. 1: Top: The order parameter k(p) for different system 
lengths L. Inset: Scaling collapse for L — 256 and larger. 
Bottom: Log-log plot of the width of the transition, W (solid 
circles) and the critical point shift 8p c (open squares) as a 
function of L with a power-law fit (line) with v = 1.52(2) for 
W and v — 1.53(1) for Sp c - The inset shows the jump k c , 
remains nonzero (solid triangles and right axis) and that p c 
remains smaller than unity (open diamonds and left axis) as 
L — > 00. 



sating force. Thus, large voids will shrink away; they 
will not grow. To capture some of this physics, we have 
introduced a 3-core model with '"force-balance" on the 
square lattice |2(| . Potential neighbors are located within 
a 5x5 square centered upon the site of interest; thus, 
Zmax = 24. To survive decimation, any given site must 
have at least 3 occupied neighbors and if there is at least 
one occupied neighbor to the right of the site of interest, 
there must be at least one occupied neighbor to its left 
and vice versa. Similarly, if there is at least one occupied 
neighbor above the site of interest, there must be at least 
one occupied neighbor below and vice versa. 

We have undertaken simulations of this 3-core model, 
and find evidence that the transition is discontinuous 
with a diverging correlation length. In Fig. la, the frac- 
tion of occupied sites in the spanning cluster, k, is plotted 
as a function of p for different system lengths L. For each 
L we observe that k jumps from zero to k c (L) at some 
p c (L). For a continuous transition, the jump k c (L) would 
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decrease with L and vanish as L — > oo, but here k c (L) 
increases with L. The size of the jump n c (L) appears to 
approach a nonzero limiting value of 0.374(1) (solid tri- 
angles in inset to Fig. lb). Furthermore, the transition 
point, p c (L) approaches p c — 0.396(1) for L — > oo. To 
verify that unstable voids do not drive the transition to 
Pc = 1, we plot 1 — p c (L) vs. 1/lni as open diamonds 
in the inset to Fig. lb. We do not find linear behavior 
with a y- intercept of zero, as predicted for p c = 1 25] . 

We calculate the correlation length exponent from two 
different quantities (Fig. lb): (1) the width of the tran- 
sition defined by W = p+{L) — p_(£), where p±(L) are 
defined as the values ofp at which the probabilities of ob- 
taining a spanning cluster are 0.25 and 0.75, respectively, 
and (2) the critical point shift Sp c = p c —p c (L). Both 
W and Sp c scale as L~ x l v with v = 1.52(2) and 1.53(1) 
respectively, as shown in Fig. lb. This exponent leads to 
scaling collapse of the order parameter curves of Fig. la, 
as shown in the inset to Fig. la. Here, we assume the 
scaling form n{p, L) = k c (L) + L~ I V f((p - p c {L))L l / v ), 
with (j = 1.0 (for optimal collapse). For ordinary first- 
order transitions, finite-size scaling would predict a di- 
verging length with an exponent of 1/d |27l ]. correspond- 
ing to 0.5 in d = 2. We obtain an very different ex- 
ponent, strongly suggesting a mixed transition. Analy- 
sis of a recently-proposed lattice model reaches a similar 
conclusion [l4|. Furthermore, a recent 1/d expansion of 
pure fc-core percolation suggests that the mixed nature 
of the mean field transition may survive in finite dimen- 
sions |2sj . 

While fc-core percolation appears to capture the mini- 
mal physics needed to explain the mixed transition found 
at Point J, it is not a complete description of jamming. 
This can already be seen by comparing the exponents ob- 
served for the mixed transition of the d = 2 3-core model, 
j3 w 1.0 and v ks 1.5, with those found in the d = 2 and 
d = 3 jamming simulations, (3 — 0.49 and v = 0.24. In 
fact, fc-core models do not include a very important prop- 
erty of the jamming transition, namely isostaticitv|l7j. 
At Point J, the number of overlapping neighbors jumps 
from zero to Z c = 2d where d is the dimensionality. In fc- 
core percolation, on the other hand, Z c is not universal; 
it depends on fc and Z max . We find that the global con- 
straint of fc-core percolation yields Z c > fc, even though 
the local constraint only requires fc neighboring occupied 
sites per site. 

Duxbury, et aL[29j have proposed that the closely- 
related problem of rigidity percolation can be mapped 
onto the fc-core percolation problem in mean-field by im- 
posing the constraint that the transition should occur 
when Z c satisfies the isostatic condition. Thus in their 
formulation, the mean-field transition occurs above p c , 
at some p r at which Z c reaches its isostatic value. They 
therefore obtain an isostatic, ordinary first-order tran- 
sition, while fc-core percolation yields a non-isostatic, 
mixed transition. Neither case properly applies to the 



jamming transition, which appears to be both mixed and 
isostatic ^2 . We note that a complementary theory 
by Wyart, et al. H3 starts with isostaticity at Point J 
and successfully describes the behavior of the density of 
vibrational modes and predicts a diverging length scale 
with exponent of 1/2, with much the same physical mean- 
ing as our 2/# = 1/2. In addition, a recent field theoret- 
ical approach also starts with isostaticity and appears 
to produce some of the same mean-field exponents that 
we find [13] • We speculate that a complete theory of 
jamming would exhibit the same mean-field exponents 
as fc- core percolation, but different behavior in finite di- 
mensions due to isostaticity; the latter effect suggests an 
upper critical dimension of two [3(1 l3ll | . 

We have argued that the physical constraint of requir- 
ing at least k = d + 1 overlapping neighbors per particle 
in zero-temperature sphere packings leads to a transition 
resembling the fc-core percolation type. However, this 
analogy may have implications ranging beyond sphere 
packings to glassforming liquids. This connection is sug- 
gested by the set of exponents we find for mean-field fc- 
core percolation, which is rare but has been seen in a 
few other models that are known to exhibit glassy dy- 
namics as the temperature is lowered. These include the 
mode-coupling theory of glasses [3411 , m ean-field theories 
of the p-component spin glass M, |3J| and kinetically- 
constrained spin models @, IS IS \Hm ■ For the latter mod- 
els, this is no coincidence since they map onto fc-core per- 
colation and its variants 0, 0, E3 • Finally, we note that 
models such as the 3-SAT spin glass are also variants of fc- 
core percolation. It can be shown that in mean field, the 
unfrustration-frustration transition has the same mixed 
character as in fc-core percolation, suggesting that the 
3-SAT spin glass may exhibit glassy dynamics. 

It has been proposed [331I36T] that the behavior of many 
jamming systems, including glasses, suspensions, foams 
and granular materials, might be captured by "jamming 
phase diagrams," in the three-dimensional space of tem- 
perature T, applied shear stress a, and packing density 
(j). In this space, the boundary separating jammed from 
unjammed behavior is a "surface" whose location is neb- 
ulous because it depends on the time-scale of the obser- 
vations, and Point J lies underneath the jamming sur- 
face. Numerical simulation results 0, E3 suggest that 
the entire jamming surface of the jamming phase dia- 
gram is indeed controlled by Point J, the unique point 
where a sharp transition occurs. Here, we have argued 
that the physics near Point J is strongly suggestive of 
the fc-core problem. The latter has a transition with un- 
usual features that mirror corresponding features found 
at Point J: a mixed first-second order transition and the 
same exponents that characterize the continuous part. In 
contrast to earlier scenarios of the glass transition based 
on avoided critical points, either at nonzero tempera- 
ture [iilj or zero temperature |H, ^ , the arguments of 
O'Hern, et al. 01 evidently suggest a scenario based on 
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an avoided mixed transition at Point J. On one hand, the 
first-order character of the transition may explain the 
presence of strong system-specific features such as the 
degree of fragility. On the other hand, the continuous 
component of the transition may explain the many ubiq- 
uitous features in the phenomenology of jamming [33j. 
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